function results=ConstructBoundaryElementK(gamma,q,coordinates)
    ls=norm(coordinates{1}-coordinates{2});
    K=zeros(length(coordinates),length(coordinates));
    for i=1:length(coordinates)
        for j=1:length(coordinates)
            if i==j
                K(i,j)=gamma*ls/6*(1+1);
            else
                K(i,j)=gamma*ls/6*(1);
            end
        end
    end
    results=K;
end